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l^"} i In this proceeding, we review three recent results. First, we show that halos formed 

f^*^ ■ in simulations with gas cooling are significantly rounder than halos formed in 

dissipationless ./V-body simulations. The increase in principle axis ratios is ~ 0.2 — 
0.4 in the inner halo and remains significant at large radii. Second, we discuss 
the CDM substructure crisis and demonstrate the sensitivity of the crisis to the 
spectrum of primordial density fluctuations on small scales. Third, we assess the 
£^ ■ ability of experiments like VERITAS and GLAST to detect 7-rays from neutralino 

' dark matter annihilation in dark subhalos about the MW. 
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1. Introduction 



A preponderance of evidence indicates that galaxies are embedded in mas- 
sive, extended dark matter (DM) halos. Simulations of structure formation 
in the hierarchical cold dark matter (CDM) paradigm predict that CDM 
halos are generally triaxial 1,2 that they teem with self-bound subhalos 3 . 

The structure of halos is an important ingredient in modeling the DM 
direct detection signals and halo shapes have recently received attention for 
testing the CDM paradigm as new and improved probes of halo shape have 
been applied 5 ' 6 . Dissipationless simulations predict that Milky Way(MW)- 
size halos have a mean minor-to-major axis ratio of c/a ~ 0.6 — 0.7 with 
a dispersion of ~ 0.1 1 , while dynamical studies suggest that the observed 
coherence of the Sagittarius tidal stream constrains the inner MW halo 
to c/a> 0.8 6 . In § 2, we present recent results on the effect of baryonic 
dissipation on halo shapes in high-resolution, cosmological simulations. 

In § 3, we turn to halo substructure. In the MW and M31, there are 
more than an order of magnitude fewer observed satellites than the pre- 
dicted number of subhalos of comparable size 3 . Several explanations have 
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Figure 1. The effect of gas cooling on halo shapes. Left: Minor-to-major axis ratio 
c/a, as a function of major axis length for a cluster-size halo. The dashed line shows 
the shape profile of the DM in the adiabatic simulation. The thick, solid line shows the 
shape profile of DM in the cooling run, while the thin, solid line shows the shape profile 
for DM and baryons in the cooling run. Right: Same as the left panel, but for a MW-size 
galaxy progenitor (see text). 



been offered, including alternative DM properties 7 and inefficient galaxy 
formation in the shallow potentials of small subhalos 8 . We study the sen- 
sitivity of the dwarf satellite population to the primordial power spectum 
(PPS) of density fluctuations on small, sub-galactic scales and demonstrate 
that our interpretation of the missing satellite problem is a function of the 
amount of small-scale power. If the lack of luminous MW satellites is due to 
inefficient galaxy formation, the MW halo should contain > 10 2 otherwise 
dark subhalos. Strong lensing will be one probe of dark subhalos 9 . More 
speculatively, the annihilation of DM particles in these dense substructures 
may result in numerous 7-ray sources in the MW halo. We assess the poten- 
tial for instruments like VERITAS 11 and GLAST 12 to detect such sources 
in favorable models of supersymmetric (SUSY) DM in § 4. 



2. Halo Shapes 

We studied the effect of gas cooling on the shapes of DM halos using high- 
resolution cosmological simulations of cluster and galaxy formation in a 
concordance ACDM cosmology. The simulations were performed with the 
ART iV-body plus Eulerian gasdynamics code 14 . We refer the reader to 
Kazantzidis et al. 13 for further details. 

Briefly, we analyzed simulations of 8 cluster-size objects of mass 
10 13 /i _1 M Q to 3 x 10 14 h^M®. The cluster simulations had a peak 
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force resolution of ~ 2.4/i _1 kpc and a DM particle mass of m p ~ 2.7 x 
10 8 /i _1 M©. We also analyzed a simulation of the early evolution (z£ 4) 
of a galaxy that becomes MW-size at z = described by Kravtsov 15 . This 
simulation had m p ~ 9.2 x 10 5 h~ 1 M Q and peak resolution ~ 183/i _1 kpc. 
The mass and force resolution are adequate to study the inner regions of 
halos reliably. For each object, we analyzed two sets of simulations started 
from the same set of initial conditions, but including different physical pro- 
cesses. In one set, the gas dynamics were treated adiabatically, without any 
radiative cooling and the results agreed well with those of iV-body simula- 
tions with no baryonic component. The second set of simulations included 
radiative cooling, and star formation. 

We measured halo shapes by diagonalizing the moment of inertia 
tensor 2 . We used "differential" shape measurements because this makes 
the axis ratios measured at each radial bin nearly independent. Our main 
results are summarized in Figure 1. In the left panel, we show the profile 
c/a, as a function of major axis length for a representative cluster-size halo. 
On the right, we show results for the galaxy progenitor. The net effect of 
baryon dissipation is striking. At small radii, the axis ratios in the cooling 
simulations are greater by A(c/a)^ 0.3 and the systematic difference per- 
sists out to ~ i? v i r , where A (c/a) ~ 0.1. The baryons in the cluster are 
mostly in a massive, central, elliptical galaxy while in the galaxy formation 
simulation, ~ 90% of the baryons are in a flattened, gaseous disk. In both 
cases the effect of cooling is weakly dependent upon radius implying that 
the effect of baryonic dissipation on halo shapes is not critically sensitive 
to the detailed morphology of the baryonic component. In addition, the 
axis ratios change with radius in a manner that is not generally monotonic, 
indicating that different regions of a system may be flattened to different 
degrees. 

3. Halo Substructure 

The most accurate technique for studying halo substructure is numeri- 
cal simulation; however, the computational expense of simulations limits 
their dynamic range and their applicability in explorations of cosmological 
parameter space. To overcome this, Zentner and Bullock (ZB) 19 devel- 
oped an approximate, analytic model for subhalo populations and an up- 
dated model has recently been successfully tested against a suite of iV-body 
simulations 20 . The model approximately accounts for the merger statistics 
of subhalos, dynamical friction, and mass loss and redistribution due to 
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Figure 2. Dwarf satellites and the power spectrum. We show the observed satellite 
velocity functions (squares) and the predicted satellite velocity functions (thick lines) 
for 6 different PS. Clockwise from the top left: standard n = 1, <rg = 0.95; n = 0.94, 
a 8 = 0.83; WMAP best-fit n = 1.03, dn/dlnfc = -0.03, ct 8 = 0.84; BSI; n = 0.84, 
ag = 0.65; and n = 0.90, erg = 0.75. The models are labeled by erg. Lines are the means 
of 100 model realizations and errorbars represent the la scatter. Observational data arc 
from the review of Mateo 21 . 

tidal forces. The model allows one to generate hundreds of realizations 
of MW-likc halos and thereby explore the distribution of possible subhalo 
populations. 

In the standard paradigm, structure forms from primordial density fluc- 
tuations characterized by a nearly scale-invariant PPS, P(k) oc k n with 
n ~ 1. This basic picture has significant observational support 16 . How- 
ever, cosmic microwave background anisotropy constrains the PPS on large 
scales, k ~ 10~ 2 foMpc -1 , while halo substructure is sensitive to small 
scale power, k ~ 10 — 100 hMpc^ 1 . ZB studied the effect of variant power 
spectra on the MW dwarf satellites. They took several PPS with various 
motivations, all normalized to COBE: (1) standard n = 1, erg = 0.95; (2) 
n = 0.94, ct§ = 0.83; (3) n = 0.9, as = 0.75; (4) running mass inflation 
n = 0.84, erg = 0.65; (5) broken scale-invariance (BSI) with a power cut-off 
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at k c = 1 /iMpc~ 117 ; and (6) the best-fit running spectrum from WMAP 
n = 1.03, dn/dlnfc = —0.03, ag = 0.84. The steps in the calculation are 
first to generate MW halo substructure realizations for each PPS and to 
model the velocity dispersions of the embedded stellar components to deter- 
mine the appropriate subhalo size (labelled by maximum circular velocity 
Vmax) in which the observed satellites may be embedded. In this way, one 
constructs predicted and observed cumulative velocity functions. 

Figure 2 summarizes the results. First, one sees that the degree to 
which the dwarf satellite problem represents a challenge is greatly alleviated 
in the WMAP best-fit cosmology. The level at which inefficient galaxy 
formation or a critical mass scale for galaxy formation must be invoked to 
solve the satellite scarcity problem is degenerate with the PPS on small 
scales. Second, the MW satellite population by itself provides independent 
evidence against extreme models, such as the low normalization, ag, = 0.65 
model which under-predict substructure. 

4. 7-rays from Dark Substructure 

One way of probing the distribution and properties of substructure as well 
as the particle nature of the DM is through the detection of gamma-rays 
from annihilations of the dark matter particle in the dense, inner regions of 
subhalos. The currently favored DM particle is provided by supersymme- 
try (SUSY) and it is the lightest of the neutralinos (x). The uncertainties 
involved in trying to deduce information about the distribution and proper- 
ties of substructure indirectly via the detection of 7-rays are twofold. First, 
there are uncertainties that stem from the underlying cosmological model 
and the details of formation of very small-scale structures 18 ' 19 and second, 
uncertainties that arise from the lack of knowledge of the mass and cou- 
plings of the dark matter particle. Using the analytic substructure model 
of § 3, we can assess the ability of experiments like VERITAS and GLAST 
to detect 7-ray fluxes from DM annihilations. 

Koushiappas et al. 22 , adopted this approach and assumed the most opti- 
mistic SUSY parameters consistent with constraints on Hm 16 to determine 
the number of expected detections at a significance S > 3, as a function 
of subhalo mass M. In order to project counts of observed subhalos be- 
yond the masses of the dwarf galaxies, several physically-motivated extrap- 
olations are necessary; however, the recent simulation of "mini-halos" at 
z ~ 26 are a first step toward justifying these extrapolations with explicit 
numerical simulations 23 . 
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Figure 3. The cumulative number of subhalos of mass M > M mln detectable at S > 3 
on the sky. Results are based on 100 realizations of a MW-size halo. Errorbars indicate 
the 68% range and down arrows indicate that > 16% realizations have zero subhalos at 
that mass. Left: The number detectable by VERITAS. The solid line shows the highest 
detection efficiency case of M x = 500GeV. For comparison, the dotted line shows results 
for M x = 200GeV and the dashed line for M x = 5TeV. Middle: The solid line shows 
our standard result for a ACDM cosmology with n = 1 and erg = 0.95. The dashed 
line shows the detectable number of subhalos with the WMAP best-fit running power 
spectrum, dn/dlnfc = —0.03. Right: The number detectable with GLAST. The dashed 
line represents the best case of M x = 50GeV in a standard ACDM cosmology. The dot- 
dashed line shows the potential number of detections for a M x = lOOGeV neutralino. 



Our results are summarized in Figure 3. The figure shows that for \ 
masses M x < lOOGeV, the large field of view of GLAST and the energy sen- 
sitivity of VERITAS will allow them to detect substructure when operated 
in concert. For example, if M x ~ 75GeV, then in the case of optimal cou- 
pling to photons there will be on average ~ 1 detectable subhalo per GLAST 
field of view. In this case, subsequent direct observations with VERITAS 
should be able to confirm the line emission feature at an energy of ~ M x 
after an exposure time of ~ 450 hours. For 100GeV< M x < 500GeV, de- 
tection requires an instrument with a large effective area, like VERITAS; 
however, such a detection must rely on serendipity due to the small num- 
ber of potentially detectable objects in VERITAS' comparably small field 
of view. For neutralino masses in excess of M x > 500GeV, substructure de- 
tection via the 7-ray signal is unlikely with either GLAST or VERITAS. 22 
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